Overcoming losses with gain in a negative refractive index metamaterial 
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On the basis of a full-vectorial three-dimensional Maxwell-Bloch approach we investigate the 
possibility of using gain to overcome losses in a negative refractive index fishnet metamaterial. We 
show that appropriate placing of optically pumped laser dyes (gain) into the metamaterial structure 
results in a frequency band where the nonbianisotropic metamaterial becomes amplifying. In that 
region both the real and the imaginary part of the effective refractive index become simultaneously 
negative and the figure of merit diverges at two distinct frequency points. 



PACS numbers: 78.67.Pt, 78.20.Ci, 42.25.Bs, 78.45.+h 

Negative refractive index metamaterials offer the pos- 
sibility of revolutionary applications, such as subwave- 
length focusing [1], invisibility cloaking [2], and "trapped 
rainbow" stopping of light [:>]. The realization of these 
materials has recently advanced from the microwave to 
the optical regime [4, 5]. However, at optical wavelengths 
metamaterials suffer from high dissipative losses due to 
the metallic nature of their constituent meta-molecules. 
It is therefore not surprising that overcoming loss restric- 
tions is currently one of the most important topics in 
metamaterials research [(»]. 

It has been suggested that, owing to causality, simul- 
taneous loss-compensation and negative refractive index 
might only be attainable in a very narrow bandwidth 
with high losses nearby [7]. In an ongoing discussion sev- 
eral authors have reasoned that causality-based criteria 
have to be applied carefully and do not in general lead to 
such a strict result [8, 9]. This said, the theoretical possi- 
bility to compensate losses in optical metamaterials does 
not necessarily imply that the gain available from opti- 
cally active media suffices to achieve this goal. Indeed, 
bulk gain coefficients are usually an order of magnitude 
smaller than the absorption coefficients of metals at op- 
tical frequencies. 

A vital clue as to how the aforementioned limita- 
tion could be overcome came from [10]. There it was 
shown that the incorporation of gain in regions of high 
field intensity gives rise to an effective gain coefficient 
that can exceed its bulk counterpart by potentially or- 
ders of magnitude. Exploiting this gain-enhancement 
effect loss reduction and giant field enhancement have 
been reported for a double-fishnet metamaterial using 
frequency-domain models [11, 12]. The particular model 
used in [11] accounts for spatial nonuniformity and sat- 
uration of the gain. However, the dynamic evolution 
of the gain system, nonlinearly pumped by a short in- 
tense pulse, cannot be described self-consistently with 
such a frequency-domain approach and relies on assump- 
tions of cw-excitation. A time-domain calculation of gain 
in two-dimensional electric or magnetic metamaterials 
has recently been reported in [13], but therein the ef- 



fect of the pump field has again not been considered self- 
consistently. 

In this Letter, we study the optical response of a neg- 
ative refractive index (NRI) metamaterial with a gain 
medium embedded in the structure and we find that com- 
plete loss-compensation and even amplification is possi- 
ble using realistic gain parameters. To this end we use a 
full-vectorial time-domain approach that manages to self- 
consistently couple the evolution of the occupation den- 
sities in the gain medium directly to Maxwell's equations 
in three dimensions [14-10]. Nonlinearity, saturation of 
the gain medium, and spatio-temporal variations of both 
absorption and emission are inherent to our model, avoid- 
ing the need for external, precalculatcd inputs. 

The considered structure is an optical double-fishnet 
metamaterial [11, 17, 18] with a square periodicity of 
p = 280 nm, perforated with rectangular holes of sides 
ttx ~ 120 nm and ay = 80 nm (see Fig. 1). The additional 
geometrical parameters are hm = 40 nm, hd = 60 nm, 
and he = 60 nm. This type of NRI metamaterial, which 




FIG. 1. (color online) . Illustration of the double-fishnet struc- 
ture with a square unit cell of side- length p highlighted. The 
two perforated silver films are embedded in a dielectric host 
material which holds the dye molecules (translucent). Dimen- 
sions are given in the text. Pump (red dashed line) and probe 
(blue solid line) pulses illustrate the pump-probe configura- 
tion with the electric field polarized along the x direction. 
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has been the topic of intense research (see, e.g., [11, 17- 
19]), exhibits low absorption compared to other meta- 
materials in the optical wavelength range. Its relatively 
low absorption makes it the most promising structure for 
complete loss-compensation [20]. 

We consider two configurations, passive and active. In 
the passive configuration two silver fishnet films are em- 
bedded inside a dielectric host that has a real refractive 
index of Uh = 1.62 (see Fig. 1). The permittivity of sil- 
ver follows a Drude model corrected by two Lorentzian 
resonances to match experimental data at visible wave- 
lengths [21]. In the active configuration we insert Rho- 
damine 800 dye molecules into the dielectric host and 
excite them optically in numerical pump-probe experi- 
ments. The chosen geometric parameters ensure a good 
overlap of the metamaterial's resonant response with the 
emission spectrum of the dye for an electric field polar- 
ization along the long side Ux of the rectangular holes 
(see Fig. 1). 

In order to self-consistently calculate the gain dynam- 
ics in this system the dye molecules are described using 
a semiclassical four-level model with two optical dipolc 
transitions [14-10]. This model is implemented by intro- 
ducing auxiliary differential equations for the position- 
and time-dependent polarization densities and oc- 
cupation densities Nj into the three-dimensional finite- 
difference time-domain (FDTD) algorithm. The time 
evolution of the polarization densities for the absorption 
{i = a) and emission (i = e) lines is then given by 
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where wo,i = (wf -f Ff)^/^ are the oscillator frequencies, 
TiLUi the electronic transition energies, F^ the half-widths 
of the resonances, and e • di the dipole strengths. The 
dye molecules embedded in the dielectric host experi- 
ence, in the Lorentz approximation, the local electric field 
Eioc = + 2) /3] E and not the average electric field 
E [22]. Saturation is accounted for by the electric field 
dependence of the occupation inversions AiVa = A^3 — iVo 
and AiVo = A^2 — -^i for absorption and emission re- 
spectively, which couple to Eq. (1). Their dynamics are 
governed by 
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the transformation from complex- to real-valued polar- 
izations ([14], p. 174). 

Absorption and emission cross-sections, taken from ex- 
perimental data, are used to calculate the dipole length 
di via (Ti = (ujQ^ie'^df /h) / (egcn/iFi), with eg being the 
vacuum permittivity and c the vacuum speed of light. 
The parameters for the four-level system are chosen as 
follows (cf. [2.H]): Ac = 27rc/wc = 710 nm, Aa = 680 nm, 
Fc = Fa l/(20fs), dc ^ 0.09nni, and da = 0.1 nm; 
T32 — Tio ~ 100 fs and T21 — 500 ps. These values 
correspond to cross-sections CTc = 2.43 x 10~^^cm^ and 
da = 3.14 X 10"-^^ cm^. We set the density of the dye 
molecules as iV = I]i=o = 6 x lO^^^cm^^ ^ 10 mM 
leading to a bulk gain coefficient of approximately g « 
N ■ (7c ^ 1460 cm" ^ at full inversion. 

In order to study the active configuration we first pump 
the dye molecules with a short, intense pulse of dura- 
tion 2ps. After a short delay of 7ps we probe the struc- 
ture with a weak broadband pulse of duration 12 fs. Fig- 
ures 2(a) and (c) show a snapshot of the spatial distribu- 
tion of the occupation inversion generated by the pump 
pulse at two perpendicular planes inside the unit cell of 
the active mctamatcrial. The effective gain coefficient 




Nonradiative decay of the occupation densities is quan- 
tified by the lifetimes Tjk- The F^Pi-terms stem from 



FIG. 2. (color online), (a) Snapshot of the occupation inver- 
sion AA'c in a plane 5nm below the upper silver fishnet film 
just before probing and (b) the electric field enhancement at 
710 nm in the same plane; both for a pump-field amplitude of 
2.0kV/cm. The silver film perforation is indicated by a white 
rectangle, (c) and (d) show the same physical quantities in a 
plane given by the propagation direction and the long side of 
the perforation. The dashed black lines in (a) and (c) high- 
light the positions of the intersection with the other plane, 
respectively. 
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in the structure can be maximized with a good matching 
between the spatial distribution of the inversion and that 
of the plasmon-enhanced electric field amplitude at the 
emission wavelength Ac [Figs. 2(b) and (d)]. Indeed, we 
see from Fig. 2 that such a matching is achieved when 
the pump and the probe have the same electric field po- 
larization. 

We remark that the considered planar structure has a 
low cavity Q- factor [Q < 50). In addition, the modal 
volume, which is directly associated with the feedback 
from metallic interfaces, is orders of magnitude larger 
than in the case of nanospheres used for "spasing" [24]. 
Therefore lasing is not expected in this structure. This is 
also supported by the fact that in our simulations we do 
not observe any appreciable depletion of the occupation 
inversion by the probe pulse. Further, on the short ps 
time-scale under investigation, which is 3 orders of mag- 
nitude shorter than the free-space spontaneous emission 
lifetime of the dye, amplified spontaneous emission is not 
expected to play a significant role [25]. 

We use the standard retrieval method [2()] to extract 
the effective refractive index n = n' + m", first, of the 
passive configuration. Note that the metamaterial struc- 
ture considered in this work is surrounded by air above 
and below the dielectric host; i.e., it is deliberately not 
placed on a thin substrate, in order to be symmetric and 
nonbianisotropic [27]. The spectral variation of n' [n") 
in the passive structure is similar to that shown by the 
cyan solid (dashed) line in Fig. 3(a) (corresponding to 
the metamaterial that has dye molecules included but is 
not pumped). We find that in this passive case the figure 
of merit FOM = —n'{X)/n"{X) has a maximum value of 
2.7 at 713 nm. 

Next, we study the active metamaterial configuration 
and, in particular, the effect that an increase of the pump 
intensity has on the obtained FOM. In the calculations 
we ensure that the gain system is probed within the 
linear regime where the standard retrieval method can 
be applied. Figure 3(a) shows the real and imaginary 
parts of the retrieved effective refractive indices for pro- 
gressively increased peak pump-field amplitudes. We see 
that increasing the pump intensity, and therefore the gain 
available for the probe, leads to a decrease in the imag- 
inary part of the refractive index in the region around 
the maximum emission cross-section of the dye (710 nm) 
indicating reduced optical losses. The stronger resonance 
that arises from the gradually intensified pump also re- 
duces the real part of the refractive index towards more 
negative values. Increasing the pump pulse amplitude 
from 0.5kV/cm to 1.5kV/cm improves and blue-shifts 
the maximum of the FOM from a value of 3 at 716 nm 
to a value of 12 at 712 nm (for the higher amplitude). 
An even higher amplitude of the pump-pulse of 2kV/cm 
results in a wavelength region (approx. 706 — 714 nm) 
where the losses in the metamaterial are completely com- 
pensated. In this region, the active metamaterial exhibits 
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FIG. 3. (color online), (a) Real and imaginary part of the re- 
trieved eflFective refractive indices of the double-fishnet struc- 
ture for different pump amplitudes. The peak electric field 
amplitude of the pump increases in steps of 0.5kV/cm from 
no pumping (cyan line, fightest) to a maximum of 2.0kV/cm 
(black line, darkest). The inset shows the real and imaginary 
part of the effective permeability (black and red line, respec- 
tively) and the result of the Kramers-Kronig relation (black 
and red dotted fines) for the highest peak electric field am- 
plitude of 2.0kV/cm. (b) The figures-of-merit (FOM) for the 
same pumping amplitudes. 



a negative absorption [cf. Figure 4(b)] and both the real 
and imaginary part of the refractive index become simul- 
taneously negative. Note from Fig. 3(b) that in this case 
the FOM diverges at the two wavelength points bound- 
ing the negative-absorption (amplification) region owing 
to n" becoming exactly zero at these two wavelengths. 

To further verify the causal nature of the obtained 
effective parameters we use the method of [28] to cal- 
culate, based on the Kramers-Kronig relations, the real 
(imaginary) part of the effective permeability from the 
imaginary (real) part of the numerically retrieved [20] 
effective permeability. An example of such a calcula- 
tion for a pump amplitude of 2 kV/cm (corresponding 
to the negative-absorption regime) is shown in the inset 
of Fig. 3(a). The excellent agreement between the results 
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FIG. 4. (color online). Detailed view of (a) the imaginary 
part of the retrieved effective refractive indices n" , (b) the 
absorption, and (c) the figures of merit (FOM) for peak pump- 
field amplitudes close to and above compensation between 1.5 
and 2.1kV/cm in steps of O.lkV/cm. 



obtained from the standard retrieval method and the 
complementary Kramers-Kronig approach further con- 
firms that the extracted parameters do obey causahty. 

Finally, Fig. 4 presents a more detailed look at n", the 
absorption coefhcient, and the FOM for pump intensi- 
ties close to and above complete loss compensation. We 
note that there is a critical amplitude of approximately 
1.85 kV/cm for the pump pulse beyond which the present 
metamaterial configuration becomes amplifying. A fur- 
ther increase of the pump field up to levels of 2.1 kV/cm 
leads to a spectral broadening of the region of negative 
absorption. At very large pump-field amplitudes above 
2.2kV/cm we observe discontinuities in the effective re- 
fractive index (data not shown here) and the Kramers- 
Kronig relation for the permeability is not obeyed any 
more. This occurs despite a smooth change in amplitude 
and phase of both the transmission and the reflection co- 
efficients with no sign of gain depletion or discontinuities 
in other physical quantities. Investigation of this regime 
will be the subject of future work. 

In conclusion, we have shown how the incorporation of 
a gain medium (Rhodamine 800 dye) into the structure of 
a double-fishnet nonbianisotropic metamaterial can fully 
compensate losses in the regime where the real part of the 



metamaterial's effective refractive index is negative. In 
this spectral range the imaginary part of the refractive 
index of the metamaterial becomes negative, too. We 
believe that this work could guide experimental efforts 
towards lossless and amplifying metamaterials, which of- 
fer access to exciting photonic applications. 
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